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I present an update on the electroweak radiative correction 
factor to semileptonic r decays, including a next-to-leading 
order resummation of large logarithms. My result differs 
both qualitatively and quantitatively from the one recently 
obtained by Davier et at. As two consequences, (i) the dis¬ 
crepancy between the predictions for the muon g — 2 based 
on r decay data and e^e“ annihilation data increases, and 
(ii) the g — 2 prediction based on r decay data appears to 
be consistent (within about one standard deviation) with the 
experimental result from BNL. 

PACS numbers: 11.10.Hi,12.15.Lk,13.35.Dx,13.40.Em. 


The largest theoretical uncertainty in the Standard 
Model prediction of the anomalous magnetic moment of 
the muon, — 2)/2, arises from the hadronic two- 

loop vacuum polarization contribution, This 

contribution is two orders of magnitude larger than the 
ultimate experimental error anticipated by the Muon 
g—2 Collaboration at BNL ||^ , so it needs to be controlled 
at the 1% level or better. While non-perturbative QCD 
effects prevent a first principles calculation, can 

be rigorously obtained experimentally from a dispersion 
relation which relates it to an integral over e“'"e“ annihi¬ 
lation cross sections. Using the conserved vector current 
(CVC) hypothesis one can obtain additional information 
by studying the invariant mass distribution of r decay 
hadronic final states. This necessitates a careful assess¬ 
ment of CVC breaking effects, which was done in a recent 
article by Davier et al. Q. In this note, I present an up¬ 
date of the short distance electroweak radiative correc¬ 
tions to T decays, representing a particular CVC breaking 
effect. This update is motivated by two mistakes in one 
of the formulas of Ref. . Numerically, the correspond¬ 
ing shifts are modest, but not negligible, and have the 
same sign. 

The leading electroweak radiative corrections to t de¬ 
cays are enhanced by a large logarithm 

'S'ew = 1 + —(1 + 2(5) In —^ = 1.01878 (1) 

where Mz = 91.1876(2) GeV [|| is the Z boson mass, and 
a = a{mr) = 1/133.50(2) Q is the QED coupling at the 
T lepton mass, rur = 1776.99(3) MeV 0, evaluated in 
the MS renormalization scheme. Q is the hypercharge of 
the weak doublet produced in the final state. Therefore, 
(5 = 1/6 for semileptonic decays, t~ h'rud{s). Since 


(5 = —1/2 for leptons, there are no large logarithms for 
leptonic r decays. 

The remaining (not logarithmically enhanced) correc¬ 
tions at 0{a) have been obtained in Ref. [|j (final state 
fermion masses are neglected throughout). In the nota¬ 
tion of Eq. (17) of Ref. they are. 
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for semileptonic and leptonic decays, respectively. In 
Ref. [Q, however, was erroneously identihed with 

the ratio. 
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This amounts to a double counting of the correction 
S'g’/^4®P _ X = —0.00432: the hadronic spectral functions 
are normalized relative to the leptonic branching ratio 
(see Eq. (10) of Ref. |^]) so that the ratio (^) must be 
included, but it is incorrect to perform an additional di- 
Vision by 5 )ew ■ 

Since numerically, 
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short distance QCD effects are of similar size as the 0{a) 
corrections discussed in the previous paragraph. They 
have been computed in Ref. and modify Eq. (p. 
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The short distance QCD correction corresponding to the 
term proportional to the strong coupling constant, a^, 
has been approximated ||] in order to obtain an analytic 
result. I have checked that this approximation repro¬ 
duces the exact 0{aas\nMl) result within about 1%. 
Since two scales enter Eq. (g), it is clear that a next-to- 
leading order renormalization group analysis is in order. 

Resummation of the leading order logarithms in 
Eq. ( 0 ) is done using the renormalization group equation 
(RGE) 
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where [3^'^ is the lowest order QED /3-function coefficient. 
This RGE is subject to the initial condition, 5'(/xo, /^o) = 
1, and its solution is given by, 







21-i?Ty 




( 8 ) 


Applied to the case at hand, this is often written as, 
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1.01937, 


where the solution to the one-loop RGE of QED, 
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has been employed. It should be stressed, that consis¬ 
tency with the RGE demands one-loop evolution of a(/r) 
within each of the factors in Eq. (!)■ On the other hand, 
the values used across the various factors, may be re¬ 
lated to each other either by one-loop evolution or includ¬ 
ing higher order running effects, since the difference is of 
higher order in the RGE (Q). The increase of S{mr, Mz) 
in Eq. i) relative to Eq. (0) due to the summation of 
©(a” In" M|) effects is about 3% of the non-resummed 
correction. 

I will now extend the RGE analysis of the previ¬ 
ous paragraph to properly sum up all logarithms of 
0(a<ln"M|). Eq. (©is' to be replaced by. 
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where /3 q ' is the lowest order QGD /3-function coefficient. 
With the definitions, 
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Eq. & is solved by, 
S{mr,Mz) = 
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= 1.01907 ±0.00001, 


where I used the solutions to the one-loop RGE of QGD, 


= -Po 


(3) a^M) 


d/i2 


(14) 


and of QEdS. The shift, -0.00030, between Eqs. (|) 
and (H) is somewhat larger than the shift, —0.00022, ob¬ 
tained in Ref. Q , which is in part due to the summation, 
but mainly due to the inputs. The uncertainty in Eq. (H) 
is from the current uncertainty in as = 0.120 ±41.002, 
while other parametric uncertainties are minisculea. Ne¬ 
glecting two-loop 0(a"+^ In" M|) effects, Eq. ( p^ ) sim¬ 
plifies. 
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which differs by only w 3 x 10“®. Neglecting further 
the numerically similar three-loop O (as a" In" M^) ef¬ 
fects, one can expand the second line in Eq. (p^ to linear 
order in a. If one then rewrites the resulting expression 
in terms of QGD scale parameters, Aqcd, one encounters 
the double-logarithmic form originally obtained 20 years 
ago 

To summarize, next-to-leading order effects reduce the 
leading order summation by about 50%, i.e., they are 
numerically of the same order. Both effects are in turn 
numerically of order a, so they must be included for a 
complete 0{a) evaluation. Unknown higher orders are 
suppressed by at least a factor of Os/tt relative to any of 
the effects mentioned before. Thus, the uncertainty due 
to higher order effects is of order O(aas) ^ 0.0003. 

Following Ref. one can write, 

r*sub,had 

5ew = S{mr,Mz) = 1-0201 ± 0.0003, (16) 

but there Sew = 1-0267 ± 0.0027 is quoted instead. Al¬ 
most 2/3 of the difference is due to the error pointed out 
after Eq. i), and about 5% due to neglecting next-to- 
leading order contributions to S{mr,Mz). Another 15% 
difference is likely due to applying Eq. (||) incorrectly (as 
discussed above). The remaining 15% can perhaps be 


^QCD corrections to Eq. ( |l0|) are suppressed by an addi¬ 
tional factor Qs/tt. Their inclusion gives rise to the summa¬ 
tion of G/ofso" In" Mz) effects, but the integration cannot be 
performed analytically. Numerically this summation affects 
the result at the 10“® level which can safely be neglected. 

^What enters Eqs. (^ and (@) is the MS 6-quark definition, 
which is free of renormalqn ambiguities and therefore much 
better known (see Ref. |ll| for a recent sub-percent determi¬ 
nation) than the 6-quarkpole mass. Higher order matching 
corrections are also smaller if one uses the MS mass definition. 


2 






































































traced to use of the on-shell definition of a in place of the 
MS definition as used in the present work. Note, that the 
derivation of Eq. (i assumes a mass-independent renor¬ 
malization scheme (such as the MS scheme), in which at 
each fermion threshold the /3-function coefficients change 
by a finite amount: this is the origin of the product 
form of Eq. (^). Thus, the solution (||) cannot be ap¬ 
plied to mass-dependent schemes, such as the on-shell 
renormalization scheme. It is emphasized again, that the 
numerical difference to Ref. Q should not be viewed as 
a scheme-dependance and thus as an estimate of uncal¬ 
culated higher order corrections (which are much smaller 
as discussed above). On the contrary, one should ex¬ 
pect that a self-consistent treatment within the on-shell 
scheme will reproduce the result of the present work. 

As far as the r-based analysis of Ref. |Q] is concerned, 
about 77% of the data is affected by S'ew- Since Sew 
obtained in this paper differs by about 0.65% from the 
one in Ref. |^, one expects a 0.5% shift in the extracted 
^^had,( 2 )^ Including an update of the CKM matrix el¬ 
ement \Vud\ entering the analysis (the value, |I4d| = 
0.9752±0.0007, is replaced by, |14d| = 0.97485±0.00046, 
from the fit result of Ref. 0), this amounts to about one 
half of the current experimental uncertainty of 0.8 parts 
per billion for the muon magnetic moment. The r- 
based Standard Model prediction would then be consis¬ 
tent with the measurement Q within about one standard 
deviation. The discrepancy to the e+e" based analysis of 
Ref. Q would correspondingly be larger. Furthermore, 
the smaller errors in Eq. (0) and in \Vud\ compared to 
Ref. [Q should lead to a slight reduction of the overall 
uncertainty of the r-based result. As a final remark, the 
recent determination [ pl| of as from the r lifetime, when 
updated with the present next-to-leading order analysis, 
increases as(Mz) by less than 0.0001. 
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